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Strange matrix elements of the nucleon 
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Results for the disconnected contributions to matrix elements of the vector current and scalar density have 
been obtained for the nucleon from the Wilson action at f3 = 6 using a stochastic estimator technique and 2000 
quenched configurations. Various methods for analysis are employed and chiral extrapolations are discussed. 



1. MOTIVATION 

Prediction, measurement and understanding of 
strangeness in the nucleon have proven to be 
significant challenges for both theory and ex- 
periment. Of particular interest at present are 
the strangeness electric and magnetic form fac- 
tors. Experimental results have been reported by 
SAMPLEjl) and HAPPEX{|, and other groups 
are also planning experiments. || 

A number of lattice QCD studies have been re- 
ported during the past few years though the 
conclusions are somewhat varied. In the present 
work we report the results of a high statistics sim- 
ulation. Different analysis techniques are stud- 
ied and chiral extrapolations are discussed. The 
strangeness scalar density and electric and mag- 
netic form factors are all analyzed together, and 
results are compared to experimental data. 

2. SIMULATIONS AND ANALYSIS 

The nucleon's strangeness form factors arise 
from a disconnected strange quark loop that we 
compute stochastically using real Zi noise Q. To 
reduce the variance, the first four terms of the 
perturbative quark matrix are subtracted. Q| (For 
the scalar density, the first five terms are sub- 
tracted.) The three matrix elements of interest 
are 



M{ S ,M,E}(t,<l) 



G 



E„ + m 
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where {N \ss\ N) = Z s G { g ] , Z s ~l- 3k 1 ,/(4k c ), 
and the Mx(t, q) can be extracted from the ratios 



R X (t,t',q) 



G^(t,t',^ 2 \t',6) 



(2) 



G( 2 Ht,0)G^(t',q) 
of 2 and 3-point correlators by various methods 
t+i 



[Rx(t,t',q) -R x (t-l,t',q)] -> Mjf(i,g),(3) 

t'=i 
t 

Y R x(t, t', q) -»■ constant + tM x (t,q), (4) 



t'=i 



Rx(t,t',q) 



constant + tMx (t, q). (5) 



Any authentic signal should be visible with each 
of these methods, and should be consistent among 
all of them.(| 

Simulations have been performed on 20 3 x 32 
lattices using the Wilson action with (3 = 6. 
Dirichlet time boundaries are used for quarks 
with the source five timesteps from the lattice 
boundary. Valence quarks use n v = 0.152, 0.153 
and 0.154. The chiral limit is k c ~ 0.1571 and our 
n v values are in the strange region. The mass of 
the quark in the loop will be held fixed at a value 
corresponding to ki = 0.152 and the analysis will 
be based on 2000 configurations with 60 noises 
per configuration. Consistent results (not shown 
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Figure 1. Ms(t,q) from Eq. (||) with K t , = m = 
0.152. From top to bottom, the plots show ql = 
0, 1, 2, 3 and 4, where = (10o/7r)g. 



in this brief article) have been obtained from the 
analysis of 100 configurations with Ki = 0.154 and 
200 Z2 noises per configuration. 

Results for M s (t,q) from Eq. (g) arc shown in 
Fig. |l} For each q considered, a clear plateau be- 
gins about 10 timesteps from the source. (The 
source is at t = in the plots.) On the other 
hand, results for Mja(t, q) and Me(i, q) are con- 
sistent with zero for each q and shown in 
Table |. 

3. CHIRAL FITS 

Quark mass and momentum dependences can 
be calculated analytically within quenched chiral 
perturbation theory from the diagrams of Fig. ^ 
as well as all tree-level counterterms. 

There are six independent parameters and we 



Table 1 

Fits to the matrix elements of Eq. @ beginning 
10 timesteps from the source, ql = (10a/ n)q. 
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0.010(82) 


0.021(81) 



consider two extreme limits: either the quenched 
77' is absent (7 = 0), or 7 ^ but the non-77' pa- 
rameters are set to zero in the loops. The physi- 
cal situation must lie between these two extremes, 
and so we expect the true values to be somewhere 
between the curves in Fig. g. (For another recent 
discussion of chiral fitting, see Ref. ||.) 

4. COMPARING TO EXPERIMENT 

The results of Figure || compare to the existing 
experimental results as follows: 



6^(^ + 0.396^) 



0.14 ±0.29 ±0.31, 
0.03 ±0.03, 

(»)/ 



Ref. 
Fig. 
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M ' 

0.025 ±0.020 ±0.014, 
0.027 ±0.016, 
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Fig. 
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where — q\ = 0.1 GeV 2 and - 
For the scalar density, we find 

(m s /m N ) (N \ss\ N) (0) = 0.15(2), 



(7) 

0.477 GeV 2 . 

(8) 



which can be compared to the result of 0.195(9) 
reported in Ref. M. 
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Figure 2. Leading loop diagrams from quenched 
chiral perturbation theory. Dashed, solid and 
double lines denote octet mesons, octet baryons 
and decuplet baryons respectively. A shaded box 
denotes a current insertion. 

We have not attempted to estimate the sizes 
of systematic uncertainties in our results due, for 
example, to quenching and to the use of chiral 
extrapolations for valence quarks in the strange 
region. The raw lattice data of Table [l] show, 
at best, only tiny strange quark effects over the 
range of momenta and quark masses that were 
studied. Large strange-quark loop contributions 
in the vector current matrix elements should be 
considered unlikely. 
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Figure 3. Lattice results extrapolated to the 
physical hadron masses. 
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